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We develop a new approach to build the eigenfunctions of a translationally shape- invariant poten- 
tial. For this we show that their logarithmic derivatives can be expressed as terminating continued 
fractions in an appropriate variable. We give explicit formulas for all the eigenstates, their specific 
form depending on the Barclay-Maxwell class to which the considered potential belongs. 
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O ' I. INTRODUCTION 

<N 

In the framework of SUSY quantum mechanics, the concept of shape invariance is a key feature of explicit exact 
solvabilityi£&i£. Combining the hierarchy of SUSY partner hamiltonians with the shape-invariance condition, it is 
(-1 ' possible, when the SUSY is unbroken, to obtain the entire spectrum of a SIP in terms of the functions characterizing 
I/"") . the partners' correspondence. Moreover it gives an access to the eigenfunctions via a generalization of the harmonic 
t-H ' creation-annihilation operators. Among all SIP those for which the partner's parameters are related by a translation 
(TSIP) play a particular role. Indeed they are so far the only ones for which we have closed-form expressions for the 
pH ' superpotentials and then for the eigenfunctions. 

In this paper we propose an alternate way to determine the eigenfunctions of a TSIP. For this we construct analyti- 
cally the logarithmic derivatives of the eigenfunctions, which we call hereafter the Riccati-Schrodinger (RS) functions 
of the system. They are solutions of a particular type of Riccati equations depending upon the energy as a parameter 
that we will call here Riccati-Schrodinger or RS equations. Using the finite-difference Backlund algorithmi^^^ 
we obtain a terminating continued fraction expression for the RS functions in terms of the superpotential (a result 
i— i, outlined in a different way and in an incomplete form by Kazimierz 11 ). Then we consider two categories of poten- 
tials which can be reduced to an harmonic or isotonic form by a change of variable satisfying a constant coefficient 
Riccati equation. We show that their ground state RS function is a first degree polynomial or a first degree Laurent 
polynomial (depending on the category under consideration) in the new variable and that they are translationally 
' shape invariant, giving general simple algebraic formula for the energy spectrum. Combining these results, we obtain 
00 . exact rational expressions for the RS functions in terms of the new variable, which permits to recover in a simple 
way the eigenfunctions of the system. Finally we establish the equivalence between these two categories and the two 
Barclay-Maxwell classes of TSI P 12 i 13 which shows that the above construction applies in fact to the whole set of TSIP. 
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A. Invariance group on the set of Riccati equations 



As established by Carinena et alii 2 -, the finite-difference Backlund algorithm is a consequence of the invariance of 

the set of Riccati equations under a subset of the group Q of smooth SL(2, M)-valued curves Map(M., SX(2,R)). For 
any element A £ Q characterized by the matrix: 



a(x) (3{x 
j(x) 

the action of A on Map(R, M) is given by 



A{x) = ( j{x) S(x) ) ' detA{x) = a{x)5{x) ~ P{x)l{x) = 1; (1) 



A a{x)w(x) + f3(x) a(x) 1 1 

w(x) -► w(x) = = —— — — rmrr - ( 2 ) 

j(x)w(x) + o(x) j(x) j(x) r y(x)w(x) + o(x) 



If A acts on a solution of the Riccati equation: 
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w'(x) = ao(x) + ai(x)w(x) + a,2(x)w (x) (3) 
we obtain a solution of a new Riccati equation: 

w'{x) — Hq{x) +ai(x)w(x) +a2(x)w 2 (x) (4) 
u 2 (x) 

the coefficients of which being given by ~u(x) = | u\(x) 

u (x) 

a (x) = M(A)~a{x) + W(x), (5) 

where: 

/ 6 2 (x) - 7 (x)S(x) 7 2 (x) \ 

M(A) = -2(3{x)6(x) a(x)S(x) + /3(x)j(x) -2a(x)j(x) , (6) 

V f3 2 {x) -a(x)(3(x) a 2 (x) J 

^ / W{i,8;x) \ 

W(x)=\ W(6,a;x) + W{(3, r ,x) (7) 
\ W(a,ftx) J 

and W(f, g;x) = f(x)g'(x) — f'(x)g(x) is the wronskian of f(x) and g(x) in x. As noted in^, Eq.© defines an affine 
action of Q on the set of general Riccati equations. 

B. Riccati-Schrodinger equations 

To a one-dimensional Schrodinger equation (ft = 1, m = |) for a potential V^a;): 

V'"(a;) + (i7-V(x))V'(x) =0 (8) 

the transformation: 

M*) - (9) 
associates a particular Riccati equation of the form: 

- w'{x) + w 2 {x) = V{x) - E (10) 

which corresponds to Eq.Q with the coefficients ao(x) = E — V(x), a\(x) — and a,2(x) = 1. 

We'll call such an equation a Riccati-Schrodinger (RS) equation and w(x) a RS function. Note that, up to a % 
factor, the RS function identifies with the quantum momentum function at energy E in the Quantum Hamilton- 
Jacobi formalism (QHJ) of Leacock and Padget t 14 i 15 From the knowledge of the RS function, we recover immediately 
the corresponding wave function via: 



ip{x) ~ exp ( — / w(s)ds ) (11) 
and the Schrodinger equation Eq([8|) can be rewritten: 
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Hip(x) = (L+L + E) i)(x) = Etp(x) 



(12) 



with 



L = - — h w(x). 
ax 



(13) 



To each node of ipix) (which is necessarily simple) corresponds a simple pole of w(x). Moreover w(x) decreases in 
the interval [xi,X2], x\ and x 2 being the turning points of the classical motion. 

C. Finite difference Backlund algorithm 

When applied to the RS equation Eq. lfTUjl . Eq© gives: 

a 2 (x) = 6 2 (x) + {E- V(x)) 7 2 (x) + W(j, 6; x) 

oi (x) = -2/3(x)5{x) - 2a(x)j(x) (E - V{x)) + W(S, a; x) + W((3, 7; x) (14) 

00 (at) = (3 2 {x) + a 2 {x) [E - V{x)) + W{a, (3; x). 

The most general elements of Q preserving the subset of RS equations has been determined in*. Among them we 
find in particular the elements of the form: 



where <f>{x) satisfies an RS equation with the same potential as Eq. fTT))) but with a shifted energy: 

-4>'{x)+tf{x) = V(x)-(E-X) 

With this choice EqfTJJ becomes 

a 2 (x) = j (-<P'{x) + 4> 2 {x) +E- V{x)) = 1 

5i(aj) = ^ + <t> 2 ( x ) + E ~ v ( x )) ~ 2( f>( x ) = 

[ a {x) = (-4>'{x) + c)> 2 (x) +E- V{x)) + X - 2(f> 2 (x) - 4>'{x) = V(x) + 2(j>'{x) 
and w(x) satisfies the RS equation: 



(15) 



(16) 



(17) 



w'(x) + w 2 (x) = V<j,{x) - X 



(18) 



where %[x) = V(x) + 2<f>'(x). 

Suppose then that the spectrum of the hamiltonian H associated to V(x) is (E^ifj^x)) with I > not necessarily 
discrete. The corresponding RS functions Wi(x) = —^j' l (x)/tp l (x) satisfy the RS equations associated with the different 
values of the energy: 



-wi(x) + wf(x) =V(x)-Ei. (19) 
Starting from a given Wk(x), for every value of / such that Ei > Ej~ , we can build a element Aki € Q of the form: 
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A k i{x) 



w k (x) Ei - E k - w k 2 (x) 



(20) 



which transforms wi as: 



I \ A ¥ ~ I \ ( \ ( \ i E l~ E k 

wi{x) -A wi{x) = w k i(x) = -w k (x) + 



w k (x) - W[(x) 



where w k i is a solution of the RS equation: 



w' kl (x)+w 2 kl {x) = V k (x)-E l 



(21) 



(22) 



with the same energy Ei as in Eq (fl3)) but with a modified potential V k (x) = V(x) + 2w' k (x) = w' k (x) + w\(x) + E k . 
This is the content of the finite-difference Backlund algorithmic^. 

In the following we use a slightly different formulation of the Backlund algorithm: 

If w k i(x) satisfy equation Eq (|2"2")) . then for every I such that Ei > E k : 



wi(x) = w k (x) 



Ei — E k 



w k (x) + w k i(x) 

In particular, choosing k = 0, we have for every excited level Ei > Eq: 



(23) 



wi(x) = w (x) - 



Ei — E a 



w (x) + w i{x) ' 



where wqi (x) is a solution of: 



with V (x) = V(x)+ 2w' (x). 



w' i(x) + wl^x) = V (x) - Ei 



(24) 



(25) 



III. RS FUNCTIONS FOR SHAPE INVARIANT POTENTIALS AS TERMINATING CONTINUED 

FRACTIONS 

A. Basics elements 

We recall some basic elements concerning SUSY quantum mechanics and shape invariance£i&2Jii. 

Let iZ_ = —d 2 /dx 2 + V-(x) be an hamiltonian the ground state energy Eo of which is supposed to be zero: Eq = 0. 
The RS function wq (x) associated to the corresponding ground state is called the superpotential of the system and 
we can write: 



F_ = L+L, 

where L = d/dx + wq(x). The supersymmetric partner of the potential is V- defined via: 



(26) 



H + = LL+ = --^+V + (x). 



(27) 



We then have: 



V±(x) = ±w' (x)+w 2 (x). 



(28) 
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The potential V-(x) is said to be shape invariant (SIP) if it depends upon a (multi)parameter a E R N and if we 
have the relation: 

V+(x,a) = V-(xJ(a))+R(a) (29) 

R (a) £ R and f(a) £ M. N being two given functions of a. 
The complete spectrum of H- is then given by: 

n-l 

E n (a) = R(a) + R( ai ) + ... + i2(o»_i) - J] R(a k ), (30) 

fc=0 

k times 

where a k — (a) = / o ... o /(a). 

When / is a simple translation f(a) = a + e, e G Mr, V- is said to be translationally shape invariant and we call 
it a TSIP. 

B. SIP and Backlund algorithm 

We suppose below that the potential of the system is an SIP, V-{x,a), with a discrete spectrum part {E n (a)) n>0 , 
the E n (a) forming an increasing sequence for every value of a with E (a) — 0. 
Eq([2S]) becomes 



- w' 0n (x,a) + wl n (x,a) = V (x,a) - E n (a) = V-(x,a) + 2w' (x,a) - E n (a), (31) 

that is 

- w' 0n (x,a) +WQ n (x,a) = V + (x,a) - E n (a). (32) 
Using the shape- invariance condition Eq ([29|) . we obtain: 



- w' 0n (x, a) + wl n (x, a) = V-(x, ai) + R(a) - E n (a), (33) 

that is, for every m 

- (w Qn (x,a) - W m (x, ax))' + wl n (x,a) - w r 2 n (x, ai) = E m (a\) ~ E n (a) + R(a). (34) 
If we take for m the specific value m = n — 1 and use Eg. ([3T)]) . Eq. (|3"4")) becomes 

- (w 0n (x,a) - Wn-iix,^))' + wl n {x, a) - w^-xix, ai) = (35) 

which is satisfied by 

w 0n (x,a) = w n -i{x,ax). (36) 

Eq ([24l) can be now rewritten as 



^R(a k ) 

w n {x,a) =w {x,a) — ^=w {x,a) — -. (37) 

wo{x,a) + w n -i(x,ai) w (x,a) + w n -i{x, ai) 
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By a direct iteration, using a standard notation for continued fractions, we arrive at: 
* For n = 1: 



w\(x, a) = wq(x, a) 



R{a) 



wq(x, a) + w Q (x, ai 



(38) 



* For n = 2: 



w 2 (x, a) = w (x, a) 



fl(o) + i?(ai) 



R(ai) 



Wq(x, a) + wq{x, ai)— wo(^, ai) + w (ie, a 2 ) 



(39) 



* For n 



R(a) + R( ai ) + gM fl(ai)+E(a 2 ) „ 

w 3 (x, a = w (x, a) — — r — — — r — — (40) 

w {x,a) + w (x,a 1 )- w (x, ai) + w (x, a 2 )- w (x, a 2 ) + w (x, a 3 ) 



and more generally: 



Wn (x, a) = M*, a) - ff\ - r . .. ^ f ... , ^.W-^-iW (41) 

tuo(x, a) + «j (x,oi)- Wo(a;, o-j-i) +w {x, a,j)— Wo(x, a n -i) + w (x, a n ) 

Inil an incomplete version of such a formula has been outlined in a very different way. Note also that, working 
to connect the QHJ formalism and SUSY quantum mechanics, Rasinariu et al.^ had obtained the recursion relation 
Eq (l37p for n = 1 in a distinct context and again in a very different way. Their exploitation of this result diverges 
from ours. 

From the knowledge of the superpotential wq(x, a), Eq. (|41[) offers, on the basis of purely algebraic manipulations, 
a direct access to the whole set of discrete excited states. It has to be compared to the known formula 7 '- giving the 
eigenstates of a SIP via the application of differential operators (generalizing the usual creator-annihilator ones) on 
the ground state: 



ip n (x, a) - L + (o)V'n-i (x, ai) ~ L + (a)...L + (a n -i)ip {x, a„) (42) 

with L(a,j) = d/dx + w (x, aj). Note that Eq l|4"2j) can be easily retrieved from the recursion relation Eq ([3T|) . Indeed 
we have from this last: 



i(x,ai) = -w (x,a) + 



E n (a) 



w Q (x,a) - w„(x,a) 



(43) 



with 



E n {a) = - {w Q (x, a) - W n (x, a))' + wl(x, a) - w^{x, oi) (44) 

Consequently 



w n -i(x,ai) = w n (x,a) 



(w (x, a) - w n (x, a))' 
wo(x,a) - w n (x,a) 



(45) 



and 



V>„_i(a:,oi) ~ (w {x,a) - w n (x,a))ip n (x,a) = L(a)if) n (x,a) (46) 

In both cases the knowledge of the ground state permits a complete reconstitution of the spectrum. Nevertheless 
Eq (f4"Tj) avoids to use successive differentiations and is invariant with respect to changes of the position variable, a 
property which is useful in the following. We see how the preceding results apply to simple systems, namely the 
one-dimensional harmonic oscillator, the isotonic oscillator, the effective radial potential for the Kepler problem and 
the Morse potential. 
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C. One dimensional harmonic oscillator 

The one dimensional harmonic oscillator potential 



V-{x) = — x 2 + V 



and the RS equation for the ground state (Eq = 0) is 



w' (x) + Wq(x) = —x 2 + V a 



(47) 



(48) 



or 



d 

dx 



(w (x) + ^x^j j (w (x) 



(49) 



Choosing Vq = — %, this equation admits clearly a polynomial solution: 



wo(x) = —x 



(50) 



which is such that w'(x) > between the turning points of the classical motion (for a bound state we must have 
in this domain E < V(x)). As required, the ground-state eigenfunction associated to this superpotential does not 
present any nodes on the real line. 
The SUSY partner of F_ (x) is then: 



V + {x)=w' Q {x)+wl{x) = — x 2 + - = V-{x)+lj 



We recognize a SIP (cf Eg ([29"])) characterized by: 



a = — = /(a), a k = a, R(a) = 2a = w 



that is a TSIP with a zero translation amplitude. 
Consequently: 



n-l 



E n (a) — Ej-i(a) — > oj = (n — j + 1) c<j 
fc=j-i 

wo(a;, a 3 _i) + wq(^, a j) = 



and Eq (|4"Tj) becomes in this case: 



iu n (x,a;) = —x 

2 wa; 



^ ^ . (n-jlllw. ,1 
[* ... f* — f* ... [* -. 



wi- 



In particular: 



(51) 



(52) 



(53) 



(54) 



w 2 (x) 



wi(x) 



LOX — 1/ X 



2ujx 



Ws(x) = %X tt-tt i~7 — \ 

°\ / 2 ujx — 2uj/{ujx — 1/x) 



3^-1) 
2 ^ i(m 2 -3) 



(55) 
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from which, with Eq (|50[) . we deduce immediately the corresponding four first eigenstates: 



ip (x) = exp (- J x wo(s)ds) ~ exp (-fx 2 ) 
ip 1 (x) = exp (— J x wi(s)ds) ~ exp (—fx 2 ) exp (f ids) = a; exp (—fa; 2 ) 

i> 2 (x) = exp (- f* w 2 (s)ds) ~ exp (-fx 2 ) exp (/* ^ds) = (^x 2 - l) exp (-fx 2 ) 



ip 3 (x) = exp (- J x w 3 (s)ds) ~ exp (-fx 2 ) exp (/" 3 ^_ 3 J ds 



3x) exp i 



(56) 



D. Isotonic oscillator 



The isotonic oscillator potential i o 17 i 18 : 



V-(x) 



J 1 2 1(1 + 1) 



4 x 2 

and the RS equation for the ground state (E = 0): 



+ V , l>0 



(57) 



w (x) +w (x) = — x + 2 



(58) 



Looking for a solution of the form wq(x) — Ax — /i/x, we obtain immediately that Eq (|58p is satisfied when A 
uj/2, (i = I + 1 and Vb = — to (I + §), giving for the superpotential: 



, , w / + 1 

w o(^) = 77^ 

2 x 

The signs have been chosen in order that w' (x) > between the turning points of the classical motion. 
The SUSY partner of V- (x) is then: 



(59) 



V+(x) = ^(x) + u)g(x) = —x 2 + i § i - -u\l + -)+u. 



We recognize a TSIP (cf Eq ([29)) ) characterized by: 



Consequently 



(f,*+l), /(o)=(|,I + 2), o fc =(|,i + l + fc), ii(a) = 2w. 



B„(a) - -Ej-iCa) = X! 2w = 2(n-j + l)w 
W (^, Oj-l) + W (x, dj) = lux — (2 (I + j) + 1) jx 



and Eq (|4"Tj) becomes in this case: 



(60) 



(61) 



(62) 



uj I + 1 2nw 
Wn(x, a )=2 X ~ — ~ ujx- (21 + 3) jx- ux-(2(l + j) + l)/x 



f f 2(n — j + f f 



2u 



ujx - (2 (I + n) - 1) jx 



(63) 



with in particular: 
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w x (x, a) = -x . (64) 

2 x lux — (21 + 3) jx 

From Eg (l59")) and Eq (|6"4")) we deduce immediately the corresponding two first eigenstates: 

^ (x) = exp (- J x w (s, a)ds) ~ exp (-fx 2 ) exp ((/ + 1) J x ±ds) = x l+1 exp (-fx 2 ) 
^(x) = exp (- J x wi (s, a)ds) ~ x l+1 exp (-fx 2 ) exp ( J x —^L^ ds) = x l+1 (cox 2 - {21 + 3)) exp i 



(65) 



E. Effective radial potential for the Kepler-Coulomb problem 

The effective radial potential for the Kepler-Coulomb i& 

V-{x) =-^ + fc^ + l/ , k>0 (66) 
x x z 

Defining y — 1/x, V— takes the form: 

V-(y) = -iv + 1(1 + l)y 2 + V , k > (67) 
and the RS equation for the ground state (Eq = 0): 

y 2 w' {y) + w 2 {y) = - 7 y + 1(1 + l)y 2 + V (68) 

Looking for a solution of the form wq(x) — — Ay + /i, we obtain immediately that Eq (|58j) is satisfied when A = 
l + l, fi = "f/2(l + 1) and Vq = 7 2 /4(Z + l) 2 , giving for the superpotential: 

«;o(i/) — (« + 1) w + 2?T+i) (69) 

The signs have been chosen in order that w'(x) > between the turning points of the classical motion. 
The SUSY partner of U_ (x) is then: 

V + (y) = ~2y 2 w' Q (y) + V- (y) = - ly + (I + 1) (I + 2) y 2 + — ^ (70) 
We recognize a TSIP (cf Eg(|2"9"))) characterized by: 



7 2 7 2 



a = 1 + 1, f(a) = 1 + 2, a k =l + l + k, R(a k ) = - — ^— . (71) 



Consequently 



4 4 4a^ +1 



n ~ 1 2 2 2 2 

£„(a) = y i?( afe ) = -JL _ JL = 2 2 (72) 

^ 4a 2 4a 2 4(/ + l) 2 4(Z + l + ?i) 2 



and 



w (x,aj-i) + w Q (x,a j ) = - (aj-i +a ] )y+ % + = -2 (l + j + ±) y + § 

Eq (|4Tj) becomes in this case: 

/ x ^ 7 7 2 /4a 2 - 7 2 /4a T 2 t 

w n (y,a) = -«V+^ - _ {a + ai)y + l/2{1/a+1/ai) _t 



... |» 



- (aj_x + aj) y + 7/2 (l/a,-_i + l/a,) - 
(an-! +a n )y + 7/2 (l/a n _i + l/a„) 



In particular: 



wi(y,a) = w (y) 



2 
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2aa\) y — j/2aai 



or 



a) = u> (a;) + 



2aai/"f — x 2aa\ 

From Eq (|69p and Eg(|76p we deduce immediately the corresponding two first eigenstates: 
ip Q (x) = exp^~y w Q (x,a)dx\ = exp (J ( - ^ ^ - — ^ dx 
~ xi+leXP ("2(I + l) a! . 



and: 



i(>i(x) = exp f — / i«x(x, a)da; I = exp I — / wq(x, a)dx I exp 



1 



x — 2aai/j 



F. Morse potential 

The Morse potential isi 

VL(x) = A 2 + B 2 e~ 2ax -2B[A+ |) e -0 *, a > 
Defining y = exp (—ax), V_ takes the form: 



V-(y) = (By-A) 2 -aBy. 
In terms of the y variable, the associated RS equation is: 
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-otyw' n {y) = E n - {By - A) 2 + aBy + w 2 n (y) 
and for the ground state, Eq = 0, we have 



ay— + (w Q {y) - (By - A)) ) (w (y) + {By - A)) = 



a solution of which is immediately obtained as 



The SUSY partner of V-(x) is then: 



w (y) = -By + A 



(82) 



(83) 



V+{y) = -ayw' {y) + u%(y) = {By - A + a) 2 - aBy + a 2 - 2a A 
We recognize a TSIP (cf Eq (|2"9")) ) characterized by: 



(84) 



= A, /(a) = A - a, R(a) = a 2 - 2a = a 2 - a\ 



Consequently 



n-l 



E n (a) = R( a k) = a — ct„ = not (2a — not) 



and 



(85) 



(86) 



n-l 



E n (a) - Ej-i(a) = ( a l ~ a l+i) = a ]-i - a l 

w {x, Oj-_i) + wo(x, dj) = a,j-i + a,j - 2By 



Eq (|4"Tj) becomes in this case: 

w n {y,a) =a- By - 
In particular: 



f ... f 



2 2 

a j-l - 0>n 



I* ... I* 



n 2 - a 1 
u n-l "n 



a + ai — 2By— a j-i + a j ~ 2_By— a n -i + a« — 2_Bj/ 



(87) 



W i(y,a) = a- J By-^ 7 ^ 



2 By 



u>2(y,a) = a- By 



a+ai— 2By— (a^— a|)/(ai +0,2— 2By) 



From Eq (f8"3"]) and Eg (155)) we deduce immediately the corresponding two first eigenstates: 



i^o(x) = exp — / wq(x, a)dx I = cxp 



a J y 



B 
a 



(89) 



and: 



ipi(x) — exp ( — / wi(x,a)dx I = exp 



a J y 



1 f Wl ^ a) dy) ~ «p ) | /J, ( 1 -n ) j (-)] i 
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IV. POTENTIALS OF FIRST AND SECOND CATEGORIES 

A. Definitions 

1. First category 

We say that a one-dimensional potential is of first category if there exists a change of variable x — > u transforming 
the potential into an harmonic one V(x) — > V(u) = X2U 2 + X\u + Ao, such that u{x) satisfies a constant coefficient 
Riccati equation: 

= A + A lU (x)+A 2 u 2 (x), (92) 

du(x)/dx being of constant sign in all the range of values of x and u. 

The one-dimensional harmonic oscillator corresponds to the special case A\ = A2 = and the Morse potential is 
generically associated to the case A\ 7^ 0, A2 — 0. These two examples have been treated above. 

We note A = A\ — AA 2 A the discriminant of the right hand side polynomial in Eg (|9"2"]l . 

If A2 7^ 0, A 7^ 0, it is always possible, via an affine change of variable u = ay + b, to reduce Eq (j9"2"|) to the canonical 
form: 

^- = a±ay 2 (x) > 0. (93) 
dx 

If A 2 7^ 0, A = 0, a straightforward affine change of variable brings Eq (|92"|) to the form: 

£=aj/ 2 (x). (94) 
that is y — -7-^- which corresponds to the radial effective Kepler-Coulomb potential studied above. 

OLyX Xq ) 

In the following we then consider that all the first category potentials are of the type: 

V(x) -> V(y) = X 2 y 2 + Xiy + Ao, A 2 > 

dy/dx — a ± ay 2 (x) > 0. ^ ' 

The new variable is then y = tan (ax + tp ) for the positive sign and y — tanh (ax + ip ) or y — coth (ax + ip a ) 
(depending on the sign of a) for the negative sign. 

In this category we find all the potentials listed in the Annexe A which contains also the three exceptional cases 
corresponding to the harmonic, Morse and effective radial Kepler-Coulomb potentials. 



2. Second category 

We say that a one dimensional potential is of second category if there exists a change of variable x — > u transforming 
the potential into an isotonic one V(x) — > V(u) = A2M 2 + Ao + JI2/U 2 , such that u(x) satisfies a constant coefficient 
Riccati equation of the form 

^ = A a + A2u 2 (x) (96) 

du(x)/dx being of constant sign in all the range of values of x and u. 

When A 2 and A are nonzero, a linear change of variable u — ay reduces Eq ([M|) to the canonical form: 

V(x) -> V(y) = \ 2 y 2 + + A , A 2 , M2 > 

dy/dx = a±ay 2 (x) > 0. ' 1 ' 

Again, the new variable is then y = tan (ax + ip ) for the positive sign and y = tanh (ax + ip ) for the negative 
sign (due to the y — ► 1/y symmetry of the functional form of V, we can limit ourselves to this two solutions). 
The example of the isotonic potential, treated above, is generically associated to the cases Ao — or A2 = 0. 
In this category we find all the potentials listed in the Annexe B. 
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V. POLYNOMIAL SOLUTIONS FOR THE RS EQUATIONS ASSOCIATED TO FIRST AND SECOND 

CATEGORY POTENTIALS 

A. First category 

The RS equation associated to a first category potential V(x) is in terms of the y variable(see Eg. ([9"S"|) ): 

2\ „,,/ 



- (a + aiy + a 2 y 2 ) w'(y) + w 2 (y) = X 2 y 2 + Xiy + A - E 
[a.Q — a — ±a 2 , cti = 0). 

We look for a polynomial solution of degree N of this equation, describing a bound state of the system. Inserting 



AT 



w(y) = Y,b n y n , b N ^0 



n=0 



into Eq. (|98)) . we obtain: 



(99) 



JV+l / N-l 2 \ 2N I N 

E E E (n+i)a m b n+1 ] y l +Y,\ E I y l = J2^- ESn ^y n - 

1=0 \ n,m=0 I n=0 
\n+m— / 



1=0 \ n=0 m=0 
\n-\-m=l 



Since ^ 0, we must have N + 1 = 2N N = 1 or 2 = 2N N = 1. The polynomial is necessarily of degree 1: 



w(y) = fro + 



(101) 



Since y(x) satisfies a constant coefficient Riccati equation, we recover the first Gendenshtein ansatz 9 . 

For a bound state E < V(x), that is w'(x) > 0, between the turning points of the classical motion and b\ is of the 
same sign as y'(x), that is positive. In this domain w(x) is not singular and the wave function does not present any 
node, which corresponds to the ground state. 

Inserting Eq. (|10ip into Eq. pOOp we obtain for the unknown coefficients bo, b\ a set of 3 equations: 



( -ah +bl = X -E 
2b Q b 1 = Ai 

I Tail + b\ = A 2 > 

which imply a constraint on E fixing the value of the ground state energy. 
After straightforward calculations we arrive at 



(102) 



b = Xi/2/3 ± (A a ) 
h = (3 ± (A a ) > 



(103) 



with 



/3 ± (A) = ±a/2 + ^{a/2) 2 + X. (104) 

The ground state energy is given by 

E = A ± aP ± (A 2 ) - (Ai/2/3 ± (A 2 )) 2 (105) 

and the corresponding RS function: 

Mv) = 7^ 1 fv-,+P±{^)y (106) 
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B. Second category 

The RS equation associated to a first category potential V(x) is in terms of the y variable(see Eq . (f97|) ) : 



(a + aiy + a 2 y 2 ) w'(y) + w 2 (y) = X 2 y 2 + A + 



y 2 



- E 



(107) 



(a = a = ±a 2 , a\ = 0). 

We look for a solution of Eq. (ll07p which is a Laurent polynomial, that is a polynomial in y and 1/y. Inserting 



N 



j (y) = E bnyn > b N, b -M + o 



--M 



into Eq. (|107p we obtain: 



(108) 



7V+1 



iV-1 



--M+1 \ n=—M-l m=0 
n+m—l 



J 



2N 

E 

l=-2M 



N 



n,m——M 
V n-\-m—l 



J 



A 2 y + A + — 
2/ 



(109) 



Since 6jv, 6_ Af ^ 0, we must have N + 1 = 2N & N = 1 or 2 = 2N & N = 1 and 2M = M - 1 M = 1 or 
2 = 2M M = 1. The Laurent polynomial is then of degree 1, as for the regular part as for the singular part: 



w(y) =b + hy + — . 

y 



(110) 



Since y{x) satisfies now a constant coefficients Riccati equation without a term of first degree, we recover the second 
Gendenshtein ansata^. 

For a bound state E < V(x), that is w'(x) > 0, between the turning points of the classical motion and b\ — b-i/y 2 
is of the same sign as y'{x), that is positive. This necessitates that &i is positive and 6_i negative. In this domain 
w(x) is not singular and the wave function does not present any node, which corresponds to the ground state. 

Inserting Eq. pT0|) into Ea.([l09]). we obtain for the unknown coefficients bo, b\ and b-% a set of 4 equations: 



' E - A = ab\ - 6q T ab-i — 2&i&_i 

bob-! = bah = 

b\ =F ha = A 2 > 

b 2 _ x + b-ia = [i 2 > 
which imply a constraint on E, fixing the value of the ground state energy. 

Eq.ljniJ gives (see Eq.lJTDJ)): 

bo = 
h = (3 ± (A 2 ) > 
6-1 = -/?+ (M 2 ) < 0. 

The ground-state energy is given by: 

Bo = A + a {(3 ± (A 2 ) ± /?+ (m 2 )) + 2/9± (A a ) /?+ (m 2 ) 
and the corresponding RS function is 



(111) 



(112) 



(113) 



(114) 
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VI. SHAPE INVARIANCE PROPERTIES OF FIRST AND SECOND CATEGORY POTENTIALS 



A. First category 

We consider a first category potential (see Eq. ([95]) ) V-(y) — A 2 y 2 +Ai?/+Ao, A2 > 0, for which dy/dx — a±ay 2 (x) > 
and E = 0. This last constraint implies (see Eq.fTUS])) that 



Ao — Aq (A2) — 



Ai 



2/3± (A 2 ) 



a0 ± (A 2 ) . 



The SUSY partner of V-{y) is 



V+(V) = V-(y) + 2w' (x) = V-(y) + M(v)^ 

= y 2 (A 2 ± 2a[3 ± (A 2 )) + X lV + A (A 2 ) + 2a(3 ± (A 2 ) . 



We make the following change of parameter: 



(115) 
(116) 



a = (3 ± (A 2 ) 



(117) 



in which case 



The initial potential is now 



and its SUSY partner is 



A-j = a (11 . fi 1 A n (" i = ( ^- I -- im. 



V-(y,a) = a (a=F a)y 2 + Aiy + A (a) 



(118) 



V+(y,a) = y 2 a (a ± a) + Xiy + Ao (a) + 2aa 

= y 2 ax (at T a) + My + A (ai) + R(a) 



(119) 



with: 



ai = a ± a 

R(a) = A (a) - A (a x ) + 2aa = (f> 1± (a) - X ± (ai) 



(120) 



where 



A 

lj± (a) = Tfl 2 + ^ 
We recognize (see Eq . (f29|) ) the characteristic behaviour of a TSIP. More generally 



(121) 



— a ± A;a 
■fl(fflfc) = </>! ± (ofc) - 1± (ajt+i) 



The spectrum (see Ea. (j3"0")0 becomes 



(122) 



-B n (a) = X ± (a) - cj> 1± (a„) 



(123) 
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E n (a) = <j> 1>± (a) ±a 2 (n±-) 2 - -*t \ (124) 
V a J 4o^ (n±a/a) 

E n is then an isotonic function of n±a/a. As for the corresponding superpotential, it takes the form (see Eq. (|106[) ): 

w (y, a) = ay + — ^, a > (125) 
Za 

Using Eq. (|4"Tj) . Eq. (|123[) and Eq. (|125p we deduce finally the general form for the RS function associated to the n th 
level of the spectrum of an arbitrary potential of the first category: 

, x , Ai <fti,± (a) - (f> 1± (q n ) 

w n (y,a) = awH ; r ■ V ... (126) 

yy ' y 2a ( + ai)tf+^(l/o + l/ai)- V ' 

<ftl, ± (gi-l) ~ <ftl,± («n) 

(oj-i + °j) y + T- (V a j-i + i/aj) - 

^ <f\± ( a n-l) - </>i,± (a n ) 

(a„ + a„_i) y + ^ (l/a n -i + V a «) 

B. Second category 

We consider a second category potential (see Eq. ([971) 1 U_(y) = A 2 y 2 + A + , for which dy/dx = a ± ay 2 (x) > 
and i^o = 0. This last constraint implies (see Eq. (|113p ) that 

A = -a (/3 ± (A 2 ) ± /3 + (/x 2 )) - 2/3 ± (A 2 ) /3 + (/i 2 ) . 

The SUSY partner of V-{y) is: 

U+(y) = U_(y) + 2 Wo (x) - V-(y) + M(v)^ (127) 
= y 2 (A 2 ± 2<*/3 ± (A 2 )) + ^ (/i 2 + 2a/3 + (/x 2 )) + A + 2a (/3 ± (A 2 ) ± /3 + (/x 2 )) 
We define the following multiparameter: 

a=(A, M ){ A = ^(A 2 ) >0, A 2 = A(A T a) 

The initial potential is now 

y-fe,a)=A(ATa)j/ 2 + A(ff : a) + A («) (129) 



with 



Its SUSY partner is 



A (o) = -a(A± n) - 2Xfj, (130) 



V+(y,a) = A (A ± a) y 2 + M (M ^ a) + A (a) + 2a (A ± /x) (131) 

2r 



V_( Vl ai)+iJ(a), 
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where 



<2i = (Xi,fii) — (X±a,fi + a) = a + A± 
R(a) = A (a) - A (01) + 2a (A ± ft) = ± (0 2 ,± (01) - ^ 2 ,± ( a )) 



(132) 



with j4± = (±a, a) and: 

2i± (a) = ^ ± (A, M ) = (A±/i) 2 - (133) 
We recognize (see Eq.([29"]l) the characteristic behaviour of a TSIP. More generally 



a,k = a + kA = (Afc,/i fc ) = (A ± ka, fi + ka) 
R(a k ) = ± (cj> 2 ± (a fe+ i) - 2 ± (afc)) , 



(134) 



The spectrum (see Ea. (|3T)|n becomes 

E n ,±(a) = ± (0 2j± (a n ) - 2 ± (a)) (135) 



or 

E n , ± (a) = ±4a 2 (n-^^j T &,± (a) (136) 

£J nj ± is then an harmonic function of n + (/i ± A) /2a. As for the corresponding superpotential, it takes the form 
(see Eq.jnUl): 

Wo{y,a) = Ay . (137) 

y 

Using Eq. (l4"Tj) , Eq. ()135[) and Eq. (|137|) , we finally deduce the general form for the RS function associated to the n th 
level of the spectrum of an arbitrary potential of the second category of the considered type: 

y (A + Ai)y- {n + H\)/yT 

f <h 2 ± ( a n) ~ 02,± (»j-l) ^ 

(Aj_i + Aj) y - + /ij) /yT 

02, ± («") ~ 02, ± («n-l) 

(A n -i + A„) y - (m„_i + Mn) /» 
It is possible to authorize a complex phase <^ in the definition of y{x) (see Eq.(|97jl) wich becomes complex. The 

real part of corresponding to a shift of the variable, we can suppose (p purely imaginary: <p = iOq. Consider the case 
where: 



d.r 
with 



a-ay 2 (x), (139) 



tanh (ax) + i tan (On) 
y(x) = tanh (ax + i9 )=-i - ° . 140 

tanh [ax) + i cot (On) 



If #o = 7r /4, l/y(x) — y(x) and the reality of iuo(y, a) and consequently of every w Hj (y, a) is ensured when fj, = — X. 
We see immediately that <f> 2 - (a), that is E n - (a) (cf Eq. (|133p and Ea. (|135p ). and V- (y, a) (cf Eq. (|129p ) are also real. 
We recover in this case the Scarf II potential (see Annexe B). 

Since tanh(# ) cannot be equal to 1, the same reasoning cannot be used when dy/dx — a + ay 2 (x), that is, when 
y{x) — tan (ax + <p ). It cannot be neither adapted for the first category potentials. 
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C. Comparison with the Barclay-Maxwell classification of TSIP 

1. Barclay-Maxwell classification of TSIP 

Concerning the case of TSIP, using exactness arguments of the SWKB condition ^ 19 ' 20 , Barclay and Maxwel l 12 ' 13 
have established that their associated supcrpotcntials obey one or other of the following equations: 
* Class I potentials: 



w' (x) = a + aiw (x) + a2wl(x) (141) 



cti 6R, i = 0,1,2 
* Class II potentials: 



w' (x) = a + aiw (x)^Ja + a 2 w^(x) + a 2 wl(x) (142) 



an 6t, i = 0,1,2 

Moreover Barclay has shown^L: 

* For Class I potentials, there exists a change of variable x — > v in which the potential is transformed into that of 
an harmonic oscillator V(x) — > V(v) = v 2 + Vo, v(x) being solution of a Riccati equation with constant coefficients: 

^^ = A + A 1 v(x)+A 2 v 2 (x) (143) 
ax 

with 



Vo = ~ , n 1 r -a . (144) 
4 (1 - a 2 ) 

* For class II potentials, there exists a change of variable x — > y in which the potential is transformed into that of 
an harmonic oscillator V(x) — > V(v) — v 2 + Vo, v(x) being solution of a first-order ODE of the form: 

dv(x) 
dx 

with 



A Q + A x v(x) sjA + A 2 v 2 (x) + A 2 v 2 (x) (145) 



T/ ao(l-a2) L a\a 2 



(146) 



2. First category 

From Eq. (|143|) and Eq. ([92|) , it appears clearly that the first category of potentials coincides with the above first 
class. Eq. (|14ip is still easily verified using the fact that the superpotential associated to a first category potential 
V(u) = X 2 u 2 + Xiu + \q is an affine function in the y variable (see Eq. (|10ip ) and then also in the u variable. 

Conversely, if w (x) obeys Eq. (|141[) . then, if we define w (x) = cu + d, c, d G M, V(x) = —w' {x) + w 2 (x) becomes 
quadratic in u and u{x) obeys a constant coefficient Riccati equation. 

3. Second category 



For the second category potential V(y) = X 2 y 2 + Ao + H 2 /y 2 , with dy(x)/dx = Aq + A 2 y 2 (x), we have obtained in 
Eq. (|137|) that the superpotential is given by: 
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w = b iy + — . (147) 

y 



Then: 



w' (x) = y'(x) ( h - = (Aoh - A 2 b. x ) + A 2 b lV 2 - A b. x \. (148) 



y J y 



From Eq. (1147|) we have also that 



V = 557 { w ° ± V w o _ 4 ^l & - 



(149) 



which combined with Eq. (|148|) gives 



with 



w' (x) = A + Bwl(x) + CJ ' A + Bwl(x) 



A = 2 (Aoh - A 2 b_ 1 ) 
B = -A/Abib-i 
I C = ± (A 6i + A 2 b- 1 ) lb x b- X ^B. 



(150) 



(151) 



Consequently any second-category potential is a Class-II potential. 

Conversely, suppose V(x) — —w' (x) + Wq(x) is a Class II potential (see Eq. (|142p ). We then define the y variable 



via: 



w = b x y ■ 



b-i 



(152) 



where &i6_i = — ao/4a 2 . Then using Eq. (|149p we have 



a 2 



= hy - 



b-i 



and from Eq. ()142j) we deduce that 



that is, 



a o , 2 



= ol 2 \ hw n + axJa 2 WQ, hw, 



Of) 



an 



2/ / V 2/ 



6-i 



(153) 



y' = A a + A 2 y 2 



(154) 



with 
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-4o = (aiy/a2~ - a 2 ) 6-1 



As for the potential, it takes the form 



V = -a - a ly fa2~wo.l— +wl + {1- a 2 )wl (156) 

V a 2 

= -a - aiVa2 (b\V + ~~ \ ~ ~~) + C 1 _ a %) f + ~^ 

y 

with 

A 2 = bf (l - a lv /a2 - a 2 ) 

M 2 = &-i(l + aiV^- a 2) (157) 
A = -ao/2 (1 + l/a 2 ) 

Consequently any Class-II potential is a second-category potential. 

Our two categories of potentials coincide then respectively with the two Barclay-Maxwell classes. Since every TSIP 

is necessarily a member of either one of these classe o 12 ' 13 , we can finally conclude that the rational representations of 
the excited states of the RS functions given in Eq. (|126[) and Eq. (|138[) are valid for the whole set of TSIP. 



(155) 



VII. CONCLUSION 



In this paper we have given an exact expression for the RS functions associated with the excited bound states of 
every SIP in terms of terminating continued fractions built on the superpotential (see Eq. (|142p b Considering the set 
of translationally shape-invariant potentials (TSIP), we have shown that it can be divided in two categories, coinciding 
with the Barclay-Maxwell classe o 12 ' 13 but using a different characterization. More precisely any TSIP can be brought 
into an harmonic or isotonic form using a change of variable resting on a constant coefficient Riccati equation without 
a term of first degree (see Eq. (p?5j) and Eq. ([W)) V It can be noticed that, according to the Chalykh-Veselov theorem, 
the harmonic and isotonic potentials are the only rational isochronous ones and they present equispaced spectr a 22 ' 23 . 
In terms of this new variable, the superpotential is a first-degree polynomial or a first-degree Laurent polynomial 
(see Eq. (|125l) and Eq. (|137p ) recovering the Gendenshtein ansatze 9 . With this formulation, the shape invariance 
characteristics of the potential appear in a very transparent manner offering a compact expression for the spectrum of 
any TSIP (see E)q. (|123[) and Eq. (|135[) ). Interestingly, we obtain a kind of duality between the TSI potentials and the 
associated energy spectrums. For a second category potential which can be brought into an isotonic form, the energy 
is an harmonic function of the shifted quantum number n ± a/a (see Eq. (|124p ). Reciprocally, for a first category 
potential which can be brought into an harmonic form, the energy is an isotonic function of n + (fi ± A) /2a (see 
Eq. (|136p ). Collecting these results together we obtain compact rational representations (see Eq. (|126p and Eq. (|138p 
for all the excited states RS functions associated to a TSIP. 
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IX. ANNEXE A 

Rosen-Morse I The Rosen-Morse I potential with zero-energy ground state is given by££: 
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Defining: 



we obtain 



with a = A and Ai = — 2B. 
The superpotential is: 



sin (ax) tan (ax) A 2 



y(x) = tan (ai — ? ) = — cot(ax), a; G ]0, ^ [ 
y'(x) = a + ay 2 (x), 



A? 



Rosen-Morse II The Rosen-Morse II potential with zero-energy ground state is given byi£: 



A A + a , , x B 2 

- — i-5 + 25 tanh (aaf) + — 7 

cosh 2 (ax) A 2 



Defining: 



y(x) = tanh (ax) , cc S 
y'(x) = a- ay 2 (x), 



we obtain: 



with a — A and Ai = 2B. 
The superpotential is: 



Eckardt The Eckardt potential with zero-energy ground state is given by 7,J 



sinh 2 (ax) A 2 

Defining: 



(158) 



V(y) = a(a - a)y 2 + \ x y + - aa, (159) 



Ai B 
wq(x) = ay + — = — Acot(ax) — — . (160) 



V(y) = ^ + 25 tanh (ax) + -j + A 2 , B < A 2 . (161) 



(162) 



^2 

V(y) = a(a + a)y 2 + X x y + - aa, (163) 



A -S 

wq(x) = ay H = Atanh(ax) H — -. (164) 

2a A 



V (g) = A{ - \ Q) - 25 coth (ax) + ^ + A 2 , B > A 2 . (165) 



y(x) = coth ((—a) x) , x G K 
y'(x) = (-a)(l-t/ 2 (x)), 



(166) 
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Poschl- Teller II The Poschl- Teller II potential with zero-energy ground state is given by^ 



V( X ) = (A-B?- MA + a) + B(B-a) 
cosh (ax) sinh (ax) 



B < A , x > 



(177) 



Defining: 



y(x) = tanh (ax) < 1 
y'(x) — a — ay 2 (x). 



(178) 



we obtain: 



V(y) =\(\ + a)y 2 + 



H(H-a) 



a(X — fx) — 2\[i, 



(179) 



with 



X = A 
li = B 



(180) 



The superpotential is: 



B 

w o(y) = Ay = A tanh (ax) - - — — — -. 

y tanh(aa;) 



(181) 



Scarf I The Scarf I potential with zero-energy ground state is given byZ&: 



Defining: 



V(x) = -A 2 + i A2 +f- A ") _ B{2A _ a) ^M 
sin (ax) sin (ax) 



, B < A, , x € 



2a 



(182) 



y(x) = tan (<f) 



y'(x) = | + f y\x). 



(183) 



we obtain: 



V^y) = A (A - a/2) j/ H 5 



-(A + A0-2A/X, 



(184) 



with A = (A + 5) /2 and fi = (A - B) /2 
The superpotential is: 



^ A + B .ax. (A-B)/2 A B 

Wq(x) = \y = taXLI ) 77-77 — = ; H ; r. 

K ' y y 2 K 2 ' tan(^f) tan(ax) sin(aa;) 



(185) 
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Scarf II The Scarf II potential with zero-energy ground state is given byZ&: 



„, , a2 (B 2 - A 2 - Aa) _ . tanh (ax) ,.,„„, 

V(x )=A 2 + ± =- — - — - — B(2A + a) -j — f. 186 

cosh 2 (ax) cosh ax 



Defining: 



y(x)=tanh(if + if ) 



!/'(!) = f - f y 2 ^). 



o j_ ii {n - a/2) _ a 
Z/ 2 2 



with A = (A + iB) /2 and ,u = (-A + iB) /2. 
The superpotential is: 



(187) 



we obtain: 

V(y) = A (A + a/2) y 2 + tl^p^L _ r ( A - M ) - 2A^, ( I s.s ) 



w (x) = \y->l = —— tanh (_+,_)_ y -———L' = A tanh M + . 189 

y 2 V 2 4/ tanh (3r + if) cosh(ax) 
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